ECTP-2011-04 



Thermodynamics of viscous Matter and Radiation in the Early universe 

A. Tawfilfl 

Egyptian Center for Theoretical Physics (ECTP), MTI University, Cairo, Egypt and 
Research Center for Einstein Physics, Freie- University Berlin, Berlin, Germany 

H. Magdy 

Egyptian Center for Theoretical Physics (ECTP), MTI University, Cairo-Egypt 

(Dated: April 20, 2012) 

Assuming that the background geometry is filled with a free gas consisting of matter and radiation 
and that no phase transitions are occurring in the early universe, we discuss the thermodynamics of 
this closed system using classical approaches. We find that essential cosmological quantities, such 
as the Hubble parameter H, scale factor a, and curvature parameter k, can be derived from this 
5— I , simple model. On one hand, it obeys the laws of thermodynamics entirely. On the other hand, 

the results are compatible with the FriedmannLemaitreRobertsonWalker model and the Einstein 
field equations. The inclusion of finite bulk viscosity coefficient derives to important changes in 
all Of these cosmological quantities. The thermodynamics of the viscous universe is studied and a 
conservation law is found. Accordingly, our picture of the evolution of the early universe and its 
astrophysical consequences seems to be the subject of radical revision. We find that the parameter 
k, for instance, strongly depends on the thermodynamics of the background matter. The time 
O ■ scale at which a negative curvature might take place, depends on the relation between the matter 

content and the total energy. Using quantum and statistical approaches, we assume that the size 
»' ! 1 of the universe is given by the volume occupied by one particle and one photon. Different types 

of interactions between matter and photon are taken into account. In this quantum treatment, 
expressions for H and a are also introduced. Therefore, the expansion of the universe turns out to 
be accessible. 
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I. INTRODUCTION 



The equations of state (EoS) describing the matter and/or radiation filling the cosmological background geometry 
play an essential role in determining the time evolution of the universe. The more reliable EoS are included, the more 
realistic the time resultant evolution is. Such models are favoured because weso farhave no observational evidence 
for the real time evolution, especially during the very early eras of the universe. The present work is an extension 
of a previous one where the cosmic background matter is conjectured to consist of two parts. The first part is a 
' massive particle with mass m. The second part is given by an absolute space (frame of reference) (i.e. that is absolute 
background) with mass M, Such a background is characterized by Newton's theory which relates the first and 
second laws to an absolute space. Rather than the mass and the gravitational field, the background itself possesses 
no other features. Seeking for completeness, we mention here that Mach reasoned that Newtons postulates, which are 
relative to the simplicity of the newtonian laws, are related to the large scale distribution of matter in the universe 
Q . Mach's postulates have been implemented in the theory of special relativity Q ■ 

In the present work, we assume that the other matter components, like dark matter and dark energy, would not 
matter much during these early stages. On the other hand, we add a new component representing radiation. The 
present treatment avoids the inclusion of the rclativistic mass of the photon. This might be the subject of a future 
work. Then, the cosmic background is to be determined by one particle with mass m and a photon with energy 
h v, where h is the Planck constant and v is the frequency. We take into consideration two cases. In the first case, 
we disregard the interactions between particle and photon. In the second case, we consider interactions, especially, 
when the quantum nature is elaborated. Like hydrogravitational dynamics, the cosmological theory Q is based on 
hypotheses similar to the ones utilized in the present paper. 

All phase transitions are disregarded and the background matter is assumed to be likely formed as a free gas (i.e. 
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nonviscous). We applied the laws of thermodynamics and the fundamentals of classical physics in order to derive 
expressions for basic cosmological quantities, such as the Hubble parameter H(t), scale factor a(t), and curvature 
parameter k. Then, we compared them with the FriedmannLcmaitrcRobcrtsonWalker (FLRW) model and Einstein's 
field equations. 

Unless it is explicitly stated, we assume natural units so that the selected universal physical constants are nor- 
malized to unity. We apply the standard cosmological model in order to gain global evidence supporting the FLRW 
model, although we disregard the relativistic and microscopic effects. The various forms of matter and radiation are 
conjectured to be homogeneously and isotropically distributed. We use non-relativistic arguments to give expressions 
for the thermodynamic quantities in the early universe, which obviously reproduce essential parts of the well-known 
FLRW model. We assume that the universe is in thermal equilibrium and therefore the interaction rates apparently 
exceed the universe expansion rate, which likely was being slowed down with an increase the comoving time t. Also, we 
assume that the expansion is adiabatic (i.e. no entropy production or heat change takes place). Finally, we take into 
consideration two forms of the cosmic background matter. The first one is an ideal gaseous fluid, which is characterized 
by the lack of interactions and the constant internal energy. The second one is a viscous fluid, which is characterized 
by the long range correlations and the velocity gradient along the scale factor a(t). 

The present paper is organized as follows. The time evolutions of the energy density in nonviscous and viscous 
cosmology are discussed in sects. [TT] and IIV1 respectively. The expansion rate itself in curved and flat universes is 
studied in section IIII1 The quantum nature of the universe is introduced in section [V] Finally, section I VII is devoted 
to discussion. 



II. RATE OF ENERGY DENSITY IN NONVISCOUS COSMOLOGY 

Based on the model introduced in Ref. we first assume that all types of energies in the early universe are heat, 
Q, and the background geometry is filled with just one particle. Then 

dQ = = dU + pdV, (1) 

where U is the internal energy and p is the pressure. The volume V can be approximated to be proportional to a 3 . It 
is obvious that ((T|) is the first law of thermodynamics. The energy density p = U/V likely decreases with the expansion 
of the universe; dp = dU /V — UdV/V 2 . In comoving coordinates, U is equivalent to the mass of the particle m. Then, 
from ([I]), we get 

dp = -3{p + p)^. (2) 

Dividing both sides by an infinitesimal time element dt results in 

p = -3(p + p)H, (3) 

which is nothing but the equation of motion from the FLRW model at a vanishing cosmological constant, A = 0, 
and curvature parameter, k = 0. The time evolution of the energy density strongly depends on the thermodynamic 
quantities, p and p, that is on the EoS of the matter and/or radiation occupying the background geometry. One dot 
means first derivative with respect to the comoving time t] H is the Hubble parameter, which relates the velocity with 
the distance, H = a/a. When inserting one photon in the background geometry, then 

P = -3(ip + p) + ^j H. (4) 

The radiation-dominated phase is usually characterized by the equation of state p = p/3. Therefore © leads to 
p (x a~ 4 . In the matter-dominated phase, p << p and therefore p cx a~ 3 , (p oc V^ 1 ). The energy density p can be 
expressed in terms of the temperature T. Then, we can rephrase the proportionality in the radiation-dominated phase 
as pocTV- 1 . 



III. EXPANSION RATE IN VISCOUS COSMOLOGY 



In previous sections, the meaning of adding a photon to the model [1| is introduced. It is supposed to represent the 
radiation filling the background geometry. Therefore, the cosmic background is now characterized by one particle with 
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mass m and one photon with frequency v. Also, it is conjectured that the photon completes one oscillation over the 
whole radius of the universe. Therefore, the photon's frequency v = 1. Such an assumption fits well with the model 
where it has been conjectured that the expansion itself is determined by the distance covered by the particle with 
mass m. In other words, the size of the universe is given by the distance covered by the particle and simultaneously 
along which the photon is able to complete one cycle. These two components are located at a distance a from some 
point in the universe. In the radial direction, the particle will have a kinetic energy ma 2 /2. The kinetic energy of 
the photon reads hv. In the opposite direction, both are affected by a gravitational force due to their masses m p , m 7 
and the mass inside the sphere, which is given as M = (47r/3)a 3 p. The latter characterizes the mass of the absolute 
background. Then, the particle's gravitational potential energy is —GMm p /a and the photon's one is —GMm^/a, 
where G is the newtonian gravitational constant. In natural units, c = ks = ft = G = 1, where c and ks are speed of 
light and the Boltzmann constant, respectively. Then, the Planck constant h = 2ir. Therefore, the total energy reads 



E 



M m r 



+ 2?r 



(5) 



where the third term represents the photo's energy. For m p 3> m-y, (O can be re- written as 

. 9 7 87T 9 

a + k = — pa . 
3 

This is nothing but the Friedman's first equation with the curvature parameter 

k = — 

m r . 



(2tt -E) 



(6) 



(7) 
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In Friedman's solution, k can be vanishing ±1, referring to flat or positively or negatively curved universe, respec- 
tively 0] ■ It is straightforward to conclude that the value assigned to k depends on the interplay between the positive 
and negative terms in (fT0| . 

For a flat universe(i.e. k = 0) it is very easy to find solutions for ([7]). One solution leads to m p — > oo (i.e. very heavy 
particle mass). The other solution relies on the photon's energy, where to the total energy is exclusively determined 
by the photon, E = 2 n. In order to omit ra p , ([7]) can be rewritten as 



2M 
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which has two solutions at k = 0, 
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(8) 



(9) 



where c is to be determined from the boundary conditions. 

It is obvious that the quantities 0(%) in ([6]) and (|7|) are coefficients and the quantity 2 it has the same dimension 
as the energy E in natural units ( ([S]) and Q, which refers to the energy of a photon that completes one cycle. For 
positively or negatively curved universe, the total energy is given by subtracting m p /2 from and adding the same 
quantity to the photons energy, respectively. In other words, the particle's mass apparently determines the curvature 
of the early universe. If it is added (to the photon's energy) it derives the universe to have closed curvature and vice 
versa. 
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(10) 



In general relatively, the curvature of space is related to the encrgymomcntum tensor. The relation is given by the 
Einstein's field equations @. 



87rT t 
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-R, 



(11) 



where R^j, and R are the Ricci curvature tensor and scalar, respectively, and g Ml/ is the metric tensor. Expression 
relates the matter (energy) content to the universe's curvature. To keep matching the assumptions of the present 
work, the cosmological constant is assumed to vanish in (jTTJ) . Equation (fTDj) . which is valid for k = ±1, leads to the 
conclusion that the curvature of the universe is positive, when the particle's mass is subtracted. When the particle's 
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mass is added, then the universe's curvature becomes negative. Obviously, such a result does not rely on a theory. 
The present work is merely designed as an effective model offering hints on the real evolution of the early universe. 

According to recent heavy-ion collision experiments @ and lattice QCD simulations matter under extreme 
conditions (very high temperature and/or pressure) seems not to be, as we used to assume over the last three decades, 
an ideal gas, in which no collisions take place. It has been found that such matter is likely fluid, that is strongly 
correlated matter with finite heat conductivity and transport properties, especially finite viscosity coefficients (bulk 
and shear) [loj]. Therefore, it is appropriate to apply this assumption on the background geometry in the early universe. 
This is the motivation of the present extension. For simplicity, we assume that the shear viscosity is almost negligible. 
Therefore, we assume that the cosmic background geometry should not necessarily be filled with an ideal free gas. In 
previous work [Tl| - [T6j . we introduced models in which we included finite bulk viscosity coefficients. The analytical 
solutions of such models are nontrivial pd| - [l6j . In the present work, we try to introduce an approach for the viscous 
cosmology using this simple model, in which we just utilize classical approaches. As will be shown below, the classical 
approach seems to work perfectly in nonviscous background matter. It is in order now to check the influences of 
viscous fluid on the cosmological evolution. The simplicity of these approaches does not enhance the validity of their 
results. Surely, it helps to come up with ideas on the reality of viscous cosmology. 

We now assume that the particle and the photon are positioned in a viscous surrounding. Then the total energy, 
([5]), gets an additional contribution from the viscosity work, which apparently would slow down the expansion of the 
universe, 

E = -m„a 2 S- - £„ a 3 - + 2n ^-f-a 3 -, (12) 

2 a a a a 

where £ is the bulk viscosity coefficient. We assume that the expansion of the universe is isotropic, that is, symmetric in 
all directions. Consequently, the shear viscosity coefficient likely vanishes. Comparing (fT12j) with Friedmann's solution 
leads to another expression for the curvature parameter, 

k = -2 2-^--a 3 H 2-^-a 3 . (13) 

m p m p a m p m p a 

For the flat universe, mp has to be very large. Otherwise, when m p remains finite, the scale factor has to be given 
as follows 

2n-E^ 1/3 



a = f 3 ^j (14) 
which apparently limits the total energy E to be less than 2 w. For positively or negatively curved universes, 



4w =F m p — 2E 
2r] 



1/3 



tV3, (15) 



respectively. In deriving (fl"4"|) and (|T5|) . the approximation a/ a ps a/ a has been applied. Otherwise, we have to find a 
solution for 

o d 2 7r — E , . 

atfc = a 3 - = , (16) 

a £ 

atfc = ±l a 3 - = 17) 

a £ 

To solve these equations, we introduce a = y. Then, a — yy a . Apparently, the solution of (jTT)) . for instance, reads 



3 2tt - E ±m p /2 



1/3 



t 1 / 3 , (18) 



2 £ 

which is almost compatible with (|15[) . When assuming that E depends on t 



2 V ? 2£ J £ J 



1/3 

(19) 
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Once again, the matter content defines the universe's curvature. The constraint on the total energy E reads, 

E > 2n T ^- (20) 

So far, we conclude that the dependence of the scale factor a on the co-moving time t significantly changes with the 
viscous property of the background matter/radiation. It is illustrated in Fig. la, where an approximate comparison 
between ideal a(t) oc t 1 / 2 and viscous background geometry a(t) oc t 1 ^ is illustrated. A comparison of the Hubble 
parameters in nonviscous and viscous background is given in Fig. lb. 




Fig. 1: (a)The dependence of the scale factor a on the co-moving time £.(b)shows the Hubble parameter H in dependence on 
t. Solid curve represents the case where the background geometry is filled with an ideal content. The dashed curve illustrates 
the effects of the viscosity. Because of the natural units, physical units are not given. 



IV. RATE OF ENERGY DENSITY IN VISCOUS COSMOLOGY 



Based on the assumptions of the present model, it seems to be allowed to include the work of the bulk viscosity, 
Eq. (P). This results in 

dU = -(pdV-(£ p +^)^dV+^ dv) . (21) 
\ a a A J 

Following the procedure given in section UH the last expression can be re-organized as done in (j6|). Then, the evolution 
of energy density, p, does not depend on the Hubble parameter H only, but also on the thermodynamic (p and p) and 
transport (£) quantities additional to, 

da 2tt 
da a 



p = -3 ( (p + p) - {, — + — ) H, (22) 



where £ = £ p + £y. Comparing this evolution equation with the one in Eckart's relativistic fluid [T7| . the relativistic 
cosmic fluid leads to a direct estimation for the bulk viscous stress n. The "conservation of total energy density" or 
its static property requires that the bulk viscous stress equals the work of bulk viscosity. 

The total thermodynamic pressure, P, is given by summing up thermodynamic and viscous pressures 

P = P+^r- (24) 



Obviously, this is one of the novel results of this model. 
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To estimate the evolution of the bulk viscous pressure, we adopt the causal evolution equation satisfying the H- 
theorem, non- negative entropy production, S l t = II 2 /£T > 0. Assuming that the total content in the background 

geometry is conserved, T?j — 0, then the rate of the energy density has to fulfil the following conservation law: 

p = -3H (Peff+P), (25) 

where p e ff = P + II. The comparison of this expression with ^ illustrates the essential effect of the bulk viscosity, 
qualitatively In order to make a quantitative comparison, we need to implement barotropic EoS, including one for EL 



V. ON THE QUANTUM COSMOLOGY 



In the previous sections, we have studied the universe as a closed system consisting of one particle, one photon 
and the absolute background. Wc have shown that this model is able to reproduce various results as the standard 
cosmological model. The quantum nature of such a system is still to be elaborated. Before doing this, some constrains 
have to be taken into account. First, we are far away from the quantization of the gravitational force [l8j |. Second, 
the gravitational constant G is taken as a universal constant, that is, it is valid always and everywhere [lj . All 
proposals and even observations about time varying G [20| are disregarded. Third, all ideas about the modification of 



the newtonian dynamics arc not implemented [21 



So far, we have checked the case of one particle and one photon in both nonviscous and viscous surroundings. 
In the following, we suggest a quantum treatment. On the one hand, it is another check for the productivity and 
projectivity of the presented model. To this destination, wc assume that the background geometry has N particles 
and photons. These quantum particles are adhered within a cubic or spherical volume, a 3 . Globally, the particles 
and photons are distributed, isotropically and homogeneously. Locally, the particles are distributed according to an 
occupation function, which depends on the particle's quantum numbers and correlations. The photons are obeying 
BoseEinstein statistics. According to the standard cosmological model, the particles and photons are allowed to 
expand in a homogeneous and isotropic way. Then, the energy of a single particle in natural unites E = (k 2 + m 2 ) 1 / 2 , 



where the momentum k = (n\x + ri2y + n^iyiix j a. To account for the interaction, we insert the potential 
correction of Uhlenbeck and Gropper [22j ■ 



J exp [1 ± exp (-mTr 2 )] dr, 



and the 



(26) 



where ± refers to fcrmions and bosons, respectively. Then, the state density in the momentum space a /h = V/(2ir) 3 , 
where a depends on t. The volume V varies with t. Based on the proposed model, the volume of the universe can be 
determined by the size that is occupied by N particles and photons, 



a(t) 3 = 4w 2 NT(t) 



g p e T <*> ^2r± 
( 

g 7 e T <*> 2r 



k 2 dk 



m p T(t)JJ ^ 



m i T(t) J J 



e T «> ± 1 
00 k 2 dk 



E-y — fl 



(27) 



where r defines the region of interaction, (i is the chemical potential, and g P (g-y) is the degeneracy factor of the particle 
(photon). The potential <j> will be introduced in section W Al 
It is apparent that differentiation with respect to the comoving time and dividing both sides by the scale factor a 
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results in the Hubble parameter 



H(t) 



where 
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7T -2rJm g T(t) g 1 
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EL 



2T(t) 



(fj, - Ej) k 2 dk 



2rJm p T{t)±^ 



9 P 



e^T(fi-E p ) u2 



gT(t) ± gT(t) 



k z dk 



- Tity 



4 r^m p T{t){4> - T{t)) ± ^(20 - 3T(t)) 



V^(3T(t) - 2c/)) + 4r^ mi T(t)(<j) - T(t)) 
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A{t) 



Pdk 



e T(t) ± e T(t 

k 2 dk 



gT(t) 



e TjtJ 



A(t) = 



2aJm p T(t) ± y/n) g p / — 
J Jo e —( 



k 2 dk 



) ± e T(t) 

k 2 dk 



T{tf 



dT(t) 
dt 



A(t) 



dT{t) 
dt 



dT(t) 
dt 

A(t) 

dT(t) 
dt 

A{t) 



(28) 



(29) 



It is essential to model dT(t)/dt (section |VB[) in order to have a numerical estimation for the time evaluation of a(t) 
and H(t). The hyperbolic trigonometric function csch(a;) = 2e x /(e 2x — 1) diverges, when x — !> 0, which is fulfilled in 
two cases; T — > oo or fi = E y . The first condition is obvious, while fi = E 7 would mean that the chemical potential 
is very small; ~ 27r, where the photon's energy and momentum arc equal. The Expression (|28|) is obviously the 
result of two contributions. The first two terms represent the contribution of the individual constituents as a free gas 
(collisionlcss). The correction of Uhlcnbcck and Gropper, ((26)) . appears in the last two terms, in which the interactions 
between the individual constituents are taken into account through 0(<p), 



A. Photoparticle interactions 



So far, we assume that the dynamics controlling the universe essentially originates in the gravitational interaction of 
the background with the particle (matter) and with the photon (radiation), separately. At ultrahigh energy, the basic 
interactions between photon and matter appear in different types. Compton scattering describes an elastic interaction. 
In natural units, the photon's energy can be given as 2w/X, where A is the wavelength. 

AE 7 = mp(l-cos(0)) -1 . (30) 

The scattering angle can be averaged as (9) « 90°. Then, the Compton scattering potential <fi cs sa m p . Summing up 
the gravitational potential with this value results in the total potential, 

4> = m p + r ^L, (31) 
a 

where m 7 is the relativistic mass of the photon. When neglecting the rclativistic mass of the photon, then Compton 
scattering results in a very small energy loss [23| , 

dE f E \ 2 _ 8tt 3 

where p 7 is the photon's density. Then, the rate of energy density reads 

P(t) = -4i*- (33) 
a D m~ 



The second type of interaction is pair production, where A + 7 — > A + e + + e~ , where the particle is participating in 
with its field. It may receive part of the energy released. The energies of a photon pair arc partly absorbed in creating 
two electrons, m e . 

2tt = 2m e + K.E, (34) 

where the lift hand side represents the photon's energy. For simplicity, it is conjectured that the total potential is 
given by the sum of the gravitational potential and the energy released from pair production, 

4> = 2U-m e + — E — -) . (35) 
V a J a 

In an astrophysical context, the energy loss through pair production [24| 

(m e T) 2 ( m e m p \ 

where the function / has been calculated in Ref. (24|. It is obvious that these two types of interactions have to be 
subtracted from the total energy, ([I2"|). (|4"3"]) and (M . 

The third type of photoparticle interactions is photopion production, where A + j^A + tt + +it~ [25|. The most 
convenient way to describe the interaction between a particle (most likely a proton, where Z = A = 1) and a photon in 
an observer's frame of reference is the invariant total energy E t in the center of momentum frame of reference (CMF), 
which moves with Lorentz factor "/cm f — E p /E t 1. The photon's energy E 1 is calculated in the proton's rest frame 
of reference [2(| . 

E t = (m 2 p + 2m p E') 1/2 , (37) 



where E = j p E 7 (l(3 p cos9) and 7 P = (1 — /^p)^ 1 / 2 . The total potential is given by 



ml + 2,n p E' 1 ) 1/2 + 2^^. (38) 



Expressions pip , (|35|1 and (|38[) can be used in (|27j) for Compton scattering, pair production, and photopion production, 
respectively. 



B. The time evolution of the temperature in the early universe 

Expression (|28[) is an essential input to estimate the energy density. For flat universe, Friedmann's solution, @, 
gives 

p = ^H{tf. (39) 
As given in (|44[) , the comoving time t is related to the energy density as follows: 

t = 6y/^p- 1/2 = 12-k^H- 1 . (40) 

Then the time derivative of temperature can be deduced as follows. 

d S = -Ull^L (41) 
dt 30F 7 dT' [ ' 

which again depends on dp/dT, the dependence of p on the cosmic temperature T. Modelling dp/dT and replacing 
(141[) in (|27|) makes it possible to prepare for an estimation for essential cosmological parameters, for instance H given 
in (|28[) . in this quantum treatment. 
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VI. DISCUSSION 



So far, we show that a classical and quantum treatment of the universe results in almost quantitatively the same 
results as the ones deduced from the standard cosmological model. The qualitative behaviour of essential cosmological 
parameters is produced. For instance, in the radiation-dominated phase, in which p = p/3. (|22[) can be solved in the 
co-moving time t. In doing this, we utilize the results obtained from the dependence of a and H on the co-moving 
time t. ThenJ5]) turns to offer a substitution of p in t. 



where q is the proportionality coefficient of a oc i 1 / 3 . This is given in (| 1 4[) or (|15[) . Expressions ([3]) and ((4]) are deduced, 
when the matter filling the background geometry is nonviscous. Under the same assumptions, the solutions of these 
two expressions, respectively, read 

PC) = i^V^ ( 43 ) 

Pit) = ^ t->. (44) 

Figure [2] shows an approximate comparison between the three cases given by (|42|) . (|43|) and (|44|) . Seeking simplicity 
we disregard all coefficients, in other words, the proportionality p oc t a is merely drawn. This is not also valid for the 
fourth term in (|42p. which has units of energy density. Therefore, no physical units can be deduced. It is obvious that 
the contribution to the energy density differs over the comoving time t. At a very early stage, the viscosity adds with 
a negligible amount to the energy density. Later on, we notice that the viscous contents seem to become dominant. 
It leads to a small increase in p with increasing t. As it is included in (|42[) with a positive sign, it apparently sets a 
limitation for the validity of the presented model. At the limit, where p increases with increasing t, the model seems 
to case being valid. 

A few comments are now in order. Along the entire history of this universe, we are assuming that the background 
geometry is filled with particles and photons. The dynamics controlling the cosmological evolution is determined by 
these two constituents, which are treated as nonviscous and viscous fluid. Furthermore, we assume that no phase 
transition is taken into consideration. Therefore, the presented model seems to assume that a certain phase remains 
unchanged, while the universe was expanding. The limitation of viscosity appearing in Fig. [2] likely would refer to the 
necessity of the phase transition accompanied by symmetry changing, for instance. 

The time evolution of the energy density has been studied in different contents filling the background geometry. 
First, we start with matter, (j3]). When adding photon (radiation), we get (j4|). The effect of the viscosity appears in 
(E2D. 

In section |Vj we make a step further towards the quantum treatment. We assume that the partition function of a 
closed system consisting of N particles and photons is able to describe the entire universe. The size of the universe 
is given by the volume occupied by these two constituents. The interactions between particle and photon are taken 
into account. All cosmological parameters, like a, H and p can be deduced. Therefore, the expansion of the universe 
seems to be accessible. 

As an outlook, we may want to check the effects of the phase transitions and the corresponding changes in the degrees 
of freedom and EoS. Also, we are planning to study the potential change when taking into account the relativistic 
mass of the photon. 
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Fig. 2: An approximate comparison between the three cases given by (|42[) and ()43|) . The solid curve represents the case where 
viscous matter and viscous radiation are filling the background geometry. The non-viscous contents are given by the dashed 
curve. Dotted curve stands for non-viscous background matter. 



[9] A. Nakamura, S. Sakai, Phys. Rev. Lett. 94, 072305 (2005); S. Sakai,, A. Nakamura, PoS LAT2005: 186, (2006); A. 

Nakamura, S. Sakai, Nucl. Phys. A 774, 775, (2006). 
[10] D. Kharzeev and K. Tuchin, JHEP 0809, 093 (2008); F. Karsch, D. Kharzeev, K. Tuchin, Phys. Lett. B 663, 217, (2008). 
[11] A. Tawfik, M. Wahba, H. Mansour and T. Harko, Annals Phys. 523, 194-207 (2011); Annalen Phys. 522, 912-923 (2010). 
[12] A. Tawfik, H. Mansour and M. Wahba, Talk given at 12th Marcel Grossmann Meeting on General Relativity, 13-18 July 

2009, Paris- France, e-Print: larXiv: 09 12.01151 [gr-qc]. 
[13] A. Tawfik, T. Harko, H. Mansour and M. Wahba, Talk at the 7th Int. Conference on Modern Problems of Nuclear Physics, 

22-25 Sep. 2009, Tashkent-Uz bekistan, Uzbek J . Phys. 12, 316-321 (2010). 
[14] A. Tawfik, T. Harko, e-Print: larXiv:1108. 56971 [astro-ph.CO] to appear in Phys. Rev. D. 
[15] A. Tawfik, Annalen Phys. 523, 423-434 (2011). 
[16] A. Tawfik and M. Wahba, Annalen Phys. 522, 849-856 (2010). 
[17] C. Eckart, Phys. Rev. 58, 919 (1940). 

[18] C. M. Will, Theory and experiment in gravitational physics: Revised edition, Cambridge University Press, Cambridge, 1993. 
[19] B.G. Sidharth, Nuovo Cim. B 115, 115 (2000). 

[20] R.K. Tiwari and U. Dwivedi, Fizika B 19, 1-8 (2010); Fizika B 19, 193-200 (2010). 
[21] R. H. Sanders, Mon. Not. Roy. Astron. Soc. 296, 1009-1018 (1998). 
[22] G. E. Uhlenbeck and L. Cropper, Phys. Rev. 41, 79 (1932). 

[23] J. L. Puget, F. W. Stecker and J. H. Bredekamp, Astrophys. J. 205, 638-654 (1976). 
[24] G. R. Blumenthal, Phys. Rev. D 1, 1596-1602 (1970). 
[25] U. Maor, Phys. Rev. 135, B1205-B1211 (1964). 

[26] A. Achterberg, Y. Gallant, C. A. Norman and D. B. Melrose, |arxiv:astro-ph/990706"o| 1999; Mannheim,K., and P. L. 
Biermann, A&A, 221, 211 (1989); V. S. Berezinskii, S. I. Grigoreva and G. T. Zatsepin, 1975, Astrophys. Space 36 3 
(1975). 



